Size of flux jumps in superconducting films 



D. V. Shantsev 1,2 , A. V. Bobyl 1,2 , Y. M. Galpcrin 1,2 , T. H. Johansen 1,3 -*, S. I. Lee 4 ' 5 
1 Department of Physics, University of Oslo, P. O. Box 1048, Blindern, 0316 Oslo, Norway 
2 A. F. Ioffe Physico-Technical Institute, Polytekhnicheskaya 26, St. Petersburg 194021, Russia 
3 Texas Center for Superconductivity and Advanced Materials, 
University of Houston, Houston, TX 77204-5002, USA 
4 National Creative Research Initiative Center for Superconductivity, Department of Physics, 
Pohang University of Science and Technology, Pohang 790-784, Korea 
5 Quantum Material laboratory, Korea Basic Science Institute, 
52 Yeoeun-Dong, Yusung-gu, Daejeon 305-333, Korea 

Magneto-optical imaging is used to visualize vortex avalanches in MgB2 films at 4 K. Avalanches 
ranging from 50 to 50000 vortices were detected. The size distribution function has a clear peak 
whose position moves towards larger sizes as the applied field increases. This field dependence as 
well as variation of flux density profile during an avalanche are well described by a proposed model 
assuming a thermal origin of the avalanches. The model is based on the adiabatic approach and 
takes into account nonlocal electrodynamics in thin superconductors. The threshold field for thermal 
avalanches is predicted to be much smaller than that for thick superconductors, in agreement with 
the experiment. 



I. INTRODUCTION 

Applying a magnetic field to a type-II superconductor 
results in formation of a metastable critical state with a 
nonuniform flux density that is sustained by vortex pin- 
ning. From an application point of view the key quan- 
tity is the critical current density that determines the 
flux density gradient. However, much more information 
about the microscopic properties can be inferred from 
the time evolution of the critical state subjected to an 
external drive, e.g. a slowly increasing applied field. The 
repulsive interaction between individual vortices driven 
through a disordered pinning landscape may then result 
in avalanche dynamics, as already observed and analyzed 
in a number of experiments, see Ref. 1 for a review. 
Some studies 2-8 have reported a power-law distribution 
of avalanche sizes that is usually interpreted as a signa- 
ture of self-organized criticality (SOC). However, in many 
experiments the distribution displays a peak at some pre- 
ferred size, 4,5,8,9 or sometimes even two peaks. 3,10,11 

These observations have motivated molecular dynam- 
ics simulations of the vortex motion, which predict a 
transition from a broad distribution to a peaked one 
as the density of pinning sites decreases. 12,13 However, 
neither the simulations nor the SOC concept take into 
account thermal effects. It is well known that vortex 
motion generates heat that makes pinning weaker and 
facilitates further motion. This positive feedback is re- 
sponsible for large catastrophic vortex avalanches, or flux 
jumps, 14 involving many millions of vortices, and usually B(x) = -B c ln 



To be able to distinguish between avalanches of ther- 
mal and non-thermal origin one needs a clear theoretical 
prediction about the size of thermal avalanches. Previ- 
ous theoretical studies of avalanche size considered slab 
superconductors placed in a parallel magnetic field. 15,16 
This geometry is dramatically different from that in most 
experiments where thin samples, often films, are placed 
in a perpendicular field. Therefore, the aim of this work 
is to 

(i) present a model which enables calculation of the size 
of flux avalanches in the perpendicular geometry, and 

(ii) compare the predictions with experimental data ob- 
tained from magneto-optical (MO) imaging of flux dy- 
namics in MgB2 films. 

We find that the avalanche size distribution, as well as 
the detailed variation of flux density profiles during an 
avalanche are in a good agreement with our model. It is 
shown that thermal flux jumps can be very small - down 
to 50 vortices - if the field is close to the threshold value. 



II. FLUX JUMP FIELD 

Consider a thin superconducting strip with edges lo- 
cated at x = ±w, thickness d, where d <C w, and infinite 
in the y direction. When the zero-field-cooled strip is 
placed in a perpendicular magnetic field B ai screening 
currents are induced the y-direction. We shall assume 
the Bean critical state, i. e., j = j c in the flux penetrated 
region, j c being the critical current density. Then the 
z-componcnt of flux density in the strip plane is ' 



x^J w 2 — a 2 + w\/x 2 — 



causes heating of the superconductor to the critical tem- 
perature. To which degree also small-scale avalanches are 
affected by thermal effects is still very unclear, and hence 
represents a major problem in understanding experimen- 
tal results. 



a < \x\ < w 



(1) 

and B(x) = in the central part of the strip, |x| < a. 
Here 



wj cosh (B a /B c ), B c = fi j c d/-!T . 



(2) 
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1: Superconducting strip in a perpendicular magnetic 



Let us assume now that a uniform fluctuation of tem- 
perature 8T occurs. Assume also that the current-voltage 
curve of the superconductor is very steep near j c , so that 
the induced flux motion proceeds faster than thermal dif- 
fusion. We can then use the adiabatic approach, i. e., as- 
sume that the heat stays where it is being released. The 
justification of this will be discussed in detail in Sec.V. 
The fluctuation leads to a decrease in the critical cur- 
rent density throughout the strip, by Sj c = \dj c /dT\ST. 
Hence, some redistribution of flux density will occur gen- 
erating the Joule heat 

X 

5Q(x) = [ jEdt=j c ( SB(x')dx', x > a. (3) 



After substituting here Eq. (1) one obtains 
djc(T) 



SQ(x) = wB c 



dT 



7(x, B a ) ST, 



with 

7{x,B a ) 
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Consider now the heat at the edge, since expectedly 
this is the most unstable region. The key quantity then is 
the ratio SQ(w)/ST, which should be compared with the 
specific heat C. If C < SQ/ST, the generated heat cannot 
be absorbed, and the temperature fluctuation will grow. 
Thus, one finds that the superconductor is unstable if 
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The solution of Eq. (6) is shown graphically in 
Fig. 2 (top). The superconductor is stable if the applied 
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FIG. 2: Stability diagrams for a strip, Eq. (6), and for a slab, 
Eq. (11). 



field is below the so-called flux jump field B^ represented 
by a solid line. For small fields, j(w, B a ) grows parabol- 
ically as 7 « B\j2B\, which allows us to find an explicit 
expression for the jump field, 



(8) 



One can estimate from Fig. 2(top) that this expression 
can be used for Bfj < B c , or, as follows from Eq. (2), 
when the flux front penetrates less than 40% of the dis- 
tance to the strip center. At high fields 7(w, B a ) in- 
creases monotonously towards the asymptotic value of 
In 2. Therefore, in a strip flux jumps never happen if 



(5 > In 2. 



(9) 



It follows from Eq. (7) that is usually a monotonously 
increasing function of temperature since C grows with T, 
while j c goes down. The temperature dependence of T* 
gives only a minor contribution possibly except for T 
very close to T c . Fig. 2 (top) can thus be considered as 
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an approximation for the stability diagram in the B-T 
plane. It means that the condition (9) defines a threshold 
temperature above which flux jumps are not observed no 
matter how large field is applied. 

For comparison, let us recall the flux jump criterion 
for a slab in a parallel field. 15,16 This problem is much 
simpler since the field profiles are linear, and the result 
is that flux jumps in a slab of width 2w can happen only 
when 
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(10) 



and the flux jump field is 



Bf h = ^w Jc ^ 



J(32 h = ^cT*, (11) 



The corresponding stability diagram is shown in 
Fig. 2 (bottom). It looks similar to that for a strip ex- 
cept for the existence of a well-defined point where the 
first-jump line ends. This difference stems from the fact 
that the applied field when the flux front reaches the mid- 
dle of the sample, the full penetration field, is finite for 
a slab, but diverges for a thin strip. 17,18 

Comparing Eq. (11) with the corresponding Eq. (8) 
for a strip one sees that the flux jump field fo r a s trip 
is smaller than that for a slab by a factor ~ */ d/w. It 
can be thought of as a "demagnetization factor" that 
characterizes the difference between the applied field B a , 
and the actual field at the strip edge. This factor is 
essentially important since the aspect ratio of the strip 
d/w is of the order of 10 for most thin-film structures. 
Hence, thin films should be much more unstable with 
respect to flux jumps than bulk superconductors. The 
threshold temperature above which flux jumps are not 
observed should also be dramatically different for films 
and for bulk samples. This is seen by comparing Eqs. (8) 
and (10) and noting that (3 is smaller than /3 s i a b by a 
factor ~ d/w. 



III. JUMP SIZE 

When the instability condition is met, any temperature 
fluctuation will trigger the development of a flux jump. 
In this section we will focus on the final state that the su- 
perconducting strip reaches after the jump. Using again 
the adiabatic approach, we find that the final temper- 
ature distribution T(x) satisfies the following equation, 
which is obtained by integration of Eq. (3), 



T{x) C{T) dT 



<Z>(x). 



(12) 



I T jc{T) 

Here To is the initial temperature before the jump, and 

<f>(x)= [ [B(x')-B (x')]dx', (13) 
Jo 



is the flux per unit length that has passed through the 
point x during the whole course of the jump. B (x) is 
the initial flux distribution. Again we employ the Bean 
critical state model, so that Bo(x) is given by Eq. (1) with 
jc = jc(To). Similarly, in the final state one has j{x) = 
j c [T(x)] for all x where B(x) ^ 0. To complete the set of 
equations, the flux and current density distributions are 
connected via the Biot-Savart law, 



B{x) = B a 



(J-od 

n 



x dx' 



(14) 



was taken into ac- 



( where the symmetry B{x) 
count). 

In their solution of the analogous case of a slab in a par- 
allel field Swartz and Bean 15 called the final state after a 
flux jump the adiabatic critical state. The adiabatic criti- 
cal state for a thin strip is therefore defined by Eqs. (12), 
(13) and (14). We shall solve these equations numerically 
for the most commonly used T dependence of j c and C, 
namely the simple linear and cubic forms, 

j c (T)=j c0 (l-T/T c ), C(T) = C(T C )(T/T C ) 3 . (15) 

Note that as the temperature increases during a jump, 
the heat capacity also grows, resulting in a stabilization 
of the flux jumps thus limiting their size. 

The solution of Eqs. (12), (13) and (14) is not neces- 
sarily unique. To ensure that we find the proper one, we 
build the solution incrementally simulating the evolution 
of the system during a flux jump. We start from a very 
small uniform temperature fluctuation ST and calculate 
the electric field E cx j(x) — j c [T(x)] everywhere where 
j > j c . The field and temperature variations are found 
as SB = dE/dx and ST = jE/C{T). Then we recalcu- 
late E, and continue until j < j c everywhere in the strip. 
The B and T distributions obtained by this procedure 
will satisfy Eqs. (12), (13) and (14) . 

Shown in Fig. 3 are flux density and temperature 
profiles in the adiabatic critical state found by solving 
Eqs. (12), (13) and (14). The flux profiles extend deeper 
into the strip than the initial critical-state profiles, which 
are shown by the dashed curves. The flux density B{x) 
inside the superconductor has mainly increased, indicat- 
ing penetration of additional flux into the strip during the 
jump. In a small region near the edge, and also on the 
outside, B(x) has decreased in order to conserve the total 
amount of flux. The temperature has increased through- 
out the region where the flux motion took place. The 
single point where B{x) did not change is obviously the 
point through which the maximum amount of flux <$>(x) 
has passed. Therefore it also shows the maximal temper- 
ature increase. 

The jump size can be characterized by the total 
amount of flux arriving into the strip during the jump, 
$(u>). Shown in Fig. 4 is <&(u>) plotted as a function of 
the applied field. For small B a there are no jumps, in full 
agreement with the instability criterion Eq. (6). As B a 
increases above the threshold, the jump size grows too, 
eventually reaching a saturation value at high fields. 
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FIG. 3: Flux density and temperature distributions before 
the jump (critical state) and after the jump (adiabatic critical 
state) for two applied magnetic fields, /3 _1 (T c /4) = In 2, and 
T = 0.1T C . 
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FIG. 4: The flux jump size as a function of applied magnetic 
field for different T and /3 _1 (T c /4) = In 2. 



IV. EXPERIMENT 

Flux jumps were studied in films of MgB2 fabricated 
on AI2O3 substrates using pulsed laser deposition. 19 The 
samples had thickness of 400 nm and lateral dimensions 
4.8x4.8 mm 2 . The films had a high degree of c-axis align- 
ment perpendicular to the plane, and showed a sharp 
superconducting transition at T c = 39 K. 

The flux density distribution was visualized using MO 
imaging based on the Faraday effect in ferrite garnet 
indicator films. For a recent review of the method, 
see Ref. 20, and a description of our setup is found 
elsewhere. 21 The sample was glued with GE varnish to 
the cold finger of the optical cryostat, and a piece of MO 
indicator covering the sample area was placed loosely on 
top of the MgB2 film. The grey levels in the MO images 
were converted to magnetic field values using a position- 
dependent calibration matrix. 

Flux penetration was studied in zero-field-coolcd films 
subjected to a slowly increasing perpendicular magnetic 
field B a . The field ramp rate was chosen sufficiently slow, 
typically 0.01 mT/s, to give results that were rate inde- 
pendent. Shown in Fig. 5(a) is a MO image of a region 
near the film edge at B a = 7.15 mT. The flux density is 
maximal at the film edge which is located near the bot- 
tom of the image. Most part of the film remains in the 
Meissner state which appears black on the image. The 
flux has here penetrated only a distance of ~ 100 /im into 
film, and the flux front is seen to be strongly nonuniform. 
Such a nonuniformity is rather typical for superconduct- 
ing films. It reflects the presence of defects at the edge 
that facilitate vortex entry. 

To detect and quantify flux jumps we shall analyze 
difference images of the type shown in Fig. 5(b-f). They 
are obtained by subtraction of subsequent MO images 
taken with 2 s intervals which implies AB a = 0.02 mT. 
The medium grey color here corresponds to unchanged 
flux density, while brighter areas show where B has in- 
creased. In all these images one can clearly see bright 
spots on the grey background, and notice that their po- 
sition always changes. These spots indicate flux jumps 
- an abrupt increase of flux density within a localized 
area. One can integrate the increase of flux density over 
the jump area to find the total amount of flux arrived 
there. E. g., for the jump shown in (g) it equals 900 
flux quanta, while the average increase of B within the 
encircled region is 2.2 mT. The duration of one jump is 
shorter than 0.1 s since the human eye could not follow 
its evolution. We emphasize that shown in Fig. 5(b-f) is 
not the full series of difference images. Approximately 
50% of images were omitted because they do not display 
any jumps. 

Shown in Fig. 6 are profiles of flux density before and 
after flux jumps. They are recorded along the x axis 
perpendicular to the film edge, as shown in Fig. 5(e). The 
curves are very similar to the theoretical ones presented 
in Fig. 3(top). Again we see an increase of B inside the 
film that is especially large close to the flux front, and a 
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FIG. 5: (a) MO image showing flux distribution near the 
edge of MgB 2 film for B a = 7.15 mT at 4 K. (b-f) Difference 
images obtained by subtracting two subsequent MO images 
taken with field interval AB a — 0.02 mT display flux jumps 
(bright spots), (g) A blow-up of the marked rectangular area 
on (f ) ; the flux arrived within the outlined jump area is 900$o , 
3?o = h/2e is the flux quantum. 



slight decrease outside. 

Let us now analyze quantitatively the size distribution 
of flux jumps. We ramped the applied field from zero to 
23 mT at 3.6 K and kept record of all flux jumps that took 
place within our observation window of 0.7 mm length 
along the film edge. The results are shown in Fig. 7 where 
every jump is represented by one symbol. The first jumps 
are detected around D„ = 3 mT and are small, i.e. con- 
taining 100 or less vortices. As B a increases, the average 
jump size steadily grows although with a large scatter of 
data. At B a = 14.5 mT we observed a jump of clearly 
different type. It is much larger than all the preceding 
jumps, extends on a distance of ~ 0.5 mm, and has a 
distinct finger-like pattern. Another jump of the same 
type was observed at B a = 17.5 mT. Similar jumps, of- 
ten called dendritic, have been reported earlier in films 
of MgB2 22-24 and several other superconductors, in par- 
ticular, YBa 2 Cu 3 07_ K , 25 ' 26 Nb, 27 Nb 3 Sn, 28 YNi 2 B 2 C, 29 
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FIG. 6: Profiles of flux density near the edge of MgB 2 film 
before and after flux jump; the jump sizes are 7500$o for 
5.6 mT and 31000<I>o for 11.6 mT.' 



Pb, 30 and NbN. 31 One sees from Fig. 7 that each den- 
dritic jump is followed by an "empty" interval of B a free 
of the small jumps. Obviously, such a massive flux in- 
trusion significantly reduces the "magnetic pressure" on 
distances of the order of dendrite length. A considerable 
increase of applied field is then needed to build up this 
pressure again in order to trigger new jumps. 

Apart from the region B a > 14 mT affected by den- 
dritic jumps, the observed & (B a ) dependence is in a good 
agreement with the curve predicted for T = 0.1T C , see 
Fig. 4. Both the experimental and theoretical Q(B a ) 
curves have a steep initial slope and tend to saturation 
at B a ~ 5Bfj. Note that this qualitative agreement is 
achieved without any fitting parameters. The only pa- 
rameter in the model, /3, is fixed by the experimental fact 
that the jumps do not occur above 10 K« T c /4. Hence, 
according to Eq. (9), we have chosen /3 _1 (T c /4) = In 2. 

The large scatter of the measured jump sizes is prob- 
ably due to sample inhomogeneities that result in a very 
nonuniform penetration depth along the flux front that 
is seen in Fig. 5(a). Another factor increasing the scatter 
is the influence of preceding avalanches that can affect 
the initial state for the following flux jumps. To quantify 
the scatter of the jump sizes we collect statistics from 
the whole sample and plot the corresponding distribu- 
tion functions in Fig. 8. The distribution strongly de- 
pends on the applied field B a , as expected from Fig. 7. 
Despite the distributions are rather broad, the increase 
of the average jump size as B a changes from 4 to 10 mT 
is quite clear. Strictly speaking, the distribution func- 
tion for B a = 4 mT might not be very accurate since we 
could have missed some jumps smaller than 50$o which 
is our resolution limit. However, for B a = 10 mT there 
is a definite peak in the range 1000-4000$o. This rules 
out any possibility for SOC behavior characterized by a 
power-law size distribution. For larger B a a distinct sec- 
ond peak appears in the size distribution due to dendritic 
avalanches. 
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FIG. 7: Statistics of flux jump sizes in MgE>2 film: every 
symbol corresponds to one jump. 




FIG. 9: ^reproducibility of flux jumps. Top: Three MO 
images taken at B a = 8.7 mT in identical experiments are 
colored as red, blue and green and merged together. Bot- 
tom: Similar color merging for difference images that show 
flux jumps that took place between B a — 8.5 and 8.7 mT. 
Positions and sizes of flux jumps are unique for every exper- 
iment, though the final flux distribution is well reproducible. 




Jump size (<J> ) 

FIG. 8: Size distribution flux jumps. The statistics is col- 
lected over all jumps that occurred in MgB2 film within 
3 mT intervals of applied field centered around B a =4, 10 and 
40 mT. 



occurred only in one experiment (green). The jump (b) 
occurred in all three experiments, however its position 
and strength were not the same. It was the strongest in 
the "blue" experiment, but somehow weaker and shifted 
to the left in the "red" experiment, and to the right in the 
"green" one. Only the jump (c) was fairly reproducible 
as seen by a purely white region. In this case the place 
of the jump was predefined by a defect in the film. The 
presence of this defect can be immediately noticed by a 
reproducible peculiar feature in the flux distribution. We 
conclude therefore that in the absence of pronounced de- 
fects a unique and irreproducible sequence of flux jumps 
in every experiment produces a well reproducible final 
flux distribution. 



V. DISCUSSION 



To study the reproducibility of the observed jumps we 
carry out several identical experiments under exactly the 
same conditions. In order to visualize the difference be- 
tween the flux distributions measured in three experi- 
ments, the three MO images were coded as red, blue and 
green, and then merged together, see Fig. 9(top). The 
image seems almost perfectly black-and-white demon- 
strating that the flux distribution is very well repro- 
ducible from one experimental run to another. However, 
if one combines in a similar way three difference images, 
the resulting image is full of colors, see Fig. 9(bottom). 
This means that positions and sizes of flux jumps that 
took place as the applied field was increasing from 8.5 to 
8.7 mT are not reproducible. For example, the jump (a) 



Most predictions of the proposed adiabatic model 
are in a good agreement with the experimental results, 
namely, (i) the field dependence of the avalanche size, 
<i>(i3 a ) shown in Figs. 4 and 7, and (ii) the flux den- 
sity profiles before and after a jump shown in Figs. 3 
and 6. It suggests that the small vortex avalanches we 
observe in the MgB2 film are of thermal origin. An- 
other argument in favor of the thermal origin is the pres- 
ence of huge dendritic jumps in the same film and in the 
same temperature range. The thermal origin of the den- 
dritic jumps is now well established both experimentally 
(they are suppressed by a contact with metal 24 ' 32 ) and 
theoretically 22,33-36 

Coexistence of rather round and elongated dendritic 
jumps, both of thermal origin, can be understood within 
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FIG. 10: Schematics of current redistribution during a uni- 
form (left) and a round (right) flux jump. The flux penetrated 
area is grey, dashed is the region where it has increased. Mag- 
netic fields generated by two infinite wires (a) and by a current 
loop (b) differ only by a factor of tt/2. 

the recent theoretical models. 33 ' 34,36 They predict three 
possible situations depending on the parameters of su- 
perconductor and the background electric field E: (i) 
stability, (ii) instability leading to uniform jumps, and 
(iii) instability leading to dendritic jumps. A transition 
from uniform to dendritic jumps should occur when E 
exceeds a threshold value. Experimentally, E is created 
by ramping magnetic field, and grows with B a . Hence, 
the appearance of dendritic jumps at larger B a seen in 
Fig. 7 perfectly agrees with the theories. 37 

An important question now is whether the proposed 
adiabatic model predicting a jump size can be applica- 
ble to jumps of any shape. The model considered an 
ideal uniform film where a jump occurs uniformly along 
the flux front. The real sample is nonuniform, and the 
instability condition cannot be simultaneously met along 
the whole length of the front. Hence, every jump remains 
localized in the y direction and eventually acquires an ap- 
proximately round shape, or grows into a long dendrite. 
We made a rough estimate of the flux density increase 
due to current redistribution corresponding to a uniform 
and to a round jump schematically shown in Fig. 10. In 
the jump center the results are different only by a numer- 
ical factor of 7r/2. The resulting temperature rise is, in 
the adiabatic approach, determined only by the local flux 
motion and, hence, independent of the jump shape. This 
justifies the use of the proposed model to predict the size 
of round jumps. However, if a jump acquires a dendritic 
shape, the current distribution changes so significantly 
that the "uniform" model fails. Strong bending of the 
current flow around the jump area creates additional ac- 
celeration of the jump development. As a result, the size 
of dendritic jumps is usually limited only by the sample 
dimensions. 

Our data show that the small jumps emerge as pre- 
cursors of dendritic jumps when the applied magnetic 
field is increasing from zero. Wc believe that this is 
a general rule and that small jumps occur everywhere 
where dendritic jumps do. This explains the presence 
of two characteristic jump sizes reported in a number of 
experiments. 3,10,11 On the other hand, one can imagine 
a situation when small flux jumps do occur, while the 
dendritic jumps do not. Then, even at very large applied 
fields there will be a special type of critical state which 



is formed by repeated local flux jumps. This was possi- 
bly the case in experiments on vortex avalanches where 
a peaked size distribution was measured. 3-5,9,10 

The adiabatic approach used in our model implies that 
magnetic flux is moving faster than heat. This is more 
than likely for the case of dendritic jumps where the 
flux is moving at velocities as high as 10-100 km/s. 25,26 
Moreover, the very fact that a jump acquires a den- 
dritic shape implies, according to recent models, 33,34,36 
that the adiabatic approach holds. At the nucleation 
stage of flux jump the flux motion must be relatively 
slower, but even here the adiabatic approach seems to 
work well. Indeed, it predicts the flux jump field, Eq. (8), 
Bfj = 1.3 mT in a reasonable agreement with experi- 
ment (for C(4K) = 0.3 kJ/Km 3 , Ref. 40). Had we been 
far from the adiabatic limit, the adiabatic model would 
have underestimated the flux jump field by a large ratio 
of thermal and magnetic diffusivities. 14 ' 39 

The total amount of flux that entered the strip during 
field increase from 10 to 11 mT can be found by subtract- 
ing two corresponding MO images. This amount turns 
out to be 2 times larger than the sum over all flux jumps 
detected in the specified range of B a . It means that only 
one half of flux arrived into the superconductor via flux 
jumps. The other half arrived via gradual penetration (or 
very small jumps that are below our resolution). This is 
not surprising. The adiabatic critical state established 
after every jump is characterized by less steep flux gradi- 
ents than the original critical state. Therefore every jump 
should be followed by a relatively quiet period when in- 
creasing B a results in a gradual penetration. The flux 
gradients build up during this period bringing the sys- 
tem to the instability threshold again. We found that 
flux jumps in the same area occur with an interval of 
AB a «0.5 mT. 11 

Our analysis suggests that thermal flux jumps in su- 
perconducting films can be virtually microscopic, down 
to at least 50 flux quanta. The thermal effects should 
therefore be carefully considered when analyzing statis- 
tics of vortex avalanches even for small avalanche sizes. 
This is especially important for Hall probe measurements 
that cannot access the spatial pattern of avalanche and 
always underestimate its actual size. 

In conclusion, we propose an adiabatic model for flux 
jumps in thin superconductors. We find the flux and 
temperature distributions in the adiabatic critical state 
after a jump, the flux jump field, and the threshold tem- 
perature above which jumps disappear. We find how 
jump size depends on applied field and demonstrate, also 
experimentally, that thermal jumps can be virtually mi- 
croscopic. In this sense they do not destroy the critical 
state, rather they present the mechanism of its formation. 
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